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ABSTRACT 

We present a study of a 20cm selected sample in the Deep SWIRE VLA Field, reaching a 5-cr limiting flux 
density at the image center of S\aghz ^ 13.5/xJy. In a 0.6 x 0.6 square degrees field, we are able to assign an 
optical/IR counterpart to 97% of the radio sources. Up to 1 1 passbands from the NUV to 4.5/im are then used to 
sample the spectral energy distribution (SED) of these counterparts in order to investigate the nature of the host 
galaxies. By means of an SED template library and stellar population synthesis models we estimate photometric 
redshifts, stellar masses, and stellar population properties, dividing the sample in three sub-classes of quiescent, 
intermediate and star-forming galaxies. We focus on the radio sample in the redshift range 0.3 < z < 1 .3 where 
we estimate to have a redshift completeness higher than 90%, and study the properties and redshift evolution of 
these sub-populations. We find that, as expected, the relative contributions of AGN and star-forming galaxies 
to the /iJy population depend on the flux density limit of the sample. At all flux levels a significant population 
of "green-valley" galaxies is observed. While the actual nature of these sources is not definitely understood, 
the results of this work may suggest that a significant fraction of faint radio sources might be composite (and 
possibly transition) objects, thus a simple "AGN vs star-forming" classification might not be appropriate to 
fully understand what faint radio populations really are. 

Subject headings: galaxies; evolution, active, starburst — radio continuum: galaxies — cosmology: observa- 
tions 



1. INTRODUCTION 

In recent years, many studies have agreed in assigning 
a relevant role to active galactic nuclei (AGN) feedback in 
shaping the evolution of galaxies, and in particular their 
star formation histories, making the co-evolution of galax- 
ies and AGNs a fundamental piece in the puzzle o f the gen- 
eral evolutio n of g a laxy populat i ons (e . g.lCroton e t al. 2006, 
Menci et al. 2006, Bowe r et al.] 120061 iMonaco et all 12007 , 
Somerville et al. 2008). As it is now believed, basically all 
massive galaxies in the local Universe harbor a massive black 
hole, and the cor relation between bl ack-hole mass and galaxy 
bulge mass (e.g. lKormendv & Gebh ardt 2001) points toward 
a close link between the formation of the black hole and of its 
host galaxy. 

At the same time, deep radio surveys have been conducted 
in association with multi-wavelength observations, allowing 
such (co-)evolution of galaxies and massive black holes to be 
probed. These deep radio surveys, for the most part at 1.4 
GHz, opened a window on the previously largely unexplored 
^Jy populations. However, unlike the Jy and mJy popula- 
tions, which are dominated by radio loud AGNs hosted by 
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quiescent galaxies, the ^Jy radio source population appears to 
be increasingly dominated by different kind of sources, star- 
forming gal axies and low-lumin osity AG Ns (e.g. among oth- 
ers Windho rst et"anil985L ICon don 19891 iJarvis & RawhnpsI 
2004). ^ 

Beside being studied at radio wa velengths (e.g. | Ciliegi et all 
1999, Richards 2000, Bon di et aLll^OOllHopkins et al]l2003l 
Huynh et aL ,2005^) , the dual nature of this composite 
fiJy populations has also been confirmed with X-r a y and 
far-infrared observations (e.g. Afonso et al.l 120011 120061 
Georgakakis et al. 2003 2004). Nonetheless, the individual 
contribution of AGNs and star-forming galaxies to the whole 
^Jy population has proved difficult to determine accurately 
for several reasons, including the often small size of the 
samples, as well as observational biases introduced, for 
instance, by optical (and in particular, but not only, spec- 
troscopic) identification of the counterparts and follow-up. 
Needless to say, this is even more true at higher redshifts, 
thus hampering our ability to set evolutionary constraints. 
Therefore, while many studies over several years have been 
devoted to this investigation, making use o f different kinds 
of information at different wavelengths (e.g. Windhorst et aU 
1985, Georgakakis et al. 1999, Gruppioni et al. 1999r 
Richards et al. 1999, Ciliegi et al. 2003, Gruppion i et al. 
| 2003. Sevmour et al. 2004, C owie et al 2004. Afonso et al. 
2005[ jHuvnh et al. 2005, Afo nso et al.l [20061 ISimpson et al. 
2006','Fomalont et al. 2006, Barger et al.'2007'.'Se vmour et al. 
2008, Ibaret al. 2008 2009, Bardelli et al. 2009)," they some- 
times have produced controversial results. 

In spite of these difficulties deep radio surveys have been 
recognized, for both the AGN and star-forming components, 
as an exceptional, powerful tool, even though their potential 
is not yet fully exploited. First, since radio emission is basi- 
cally unaffected by dust extinction, important issues at opti- 
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cal wavelengths, e.g. obscured star formation and highly ob- 
scured AGNs r nissing from deep X-ray surveys (b ut see dis- 
cussions m e.g. iBarger et alJlIOOTllTozzi et alj |2009, and ref- 
erences therein), clearly find a solution when observing at ra- 
dio wavelengths. In fact, radio-selected AGN samples are not 
the same as AGN samples selected at other wavelengths, since 
they include populations of low-power radio sources which 
would not be class ified as AGNs from their optical or X- 
ray properties (e.g . iBest etalj|2005al iHardcastle etal.1 120061 
iHickox et al.ll200 91. pointing toward an intrinsically different 
nature of these sources. Furthermore, the arcsecond resolu- 
tion available for some of these radio surveys makes it rela- 
tively easy to cross-correlate them with other data across a 
broad wavelength range including optical and near-infrared. 
This is actually a fundamental point, because in fact the study 
of the faint radio sources at many different wavelengths ob- 
viously maximizes the scientific return of the radio survey, 
allowing a more complete characterization of a population 
which is intrinsically mixed at radio wavelengths. In partic- 
ular, the cross-correlation with large X-ray/optical/NIR sur- 
veys where a wealth of information is available in terms of 
spectroscopic/photometric redshifts, stellar populations and 
galaxy morphologies, enhances our understanding of the na- 
ture of these sources, out to redshift « 1. Just in the last 
couple of years, several studies were published making use 
of such deep, panchromatic observations in order to inves- 
tigate the different galaxy species populating the fjJy sam- 
ples, as for inst ance Smolcic et al. (2008) in t he CO SMOS 
field , Mainieri et al.1 (12008) and Padovani et ak (2009) in the 
GOODS-CDFS field, and Huynh et al. (2008) in the HDF-S. 

This paper is the fourth in a series documenting our 
study of the deep SWIRE field centere d at 10H6" '00', 
59°01'00"(J2000). Paper I ( Owen & MoSio^ |200i) de- 
scribes the 20cm VLA observations which produced the deep- 
est 20cm radio survey to date with 2050 sources and the 
basic radio prop erties of the faint /iJy population. Paper II 
(lOwen et alJl2009) details a complementary, deep 90cm sur- 
vey and dependence of 20cm to 90cm spectral ind ex on radio 
flux density. Paper III ( Owen & Morrisoiil 120091) documents 
the WIYN spectroscopy of sources in this field. This paper 
deals with the first properties derived for the ^Jy population, 
namely the photometric redshifts and inferred redshift distri- 
bution, and the stellar population properties of the host galax- 
ies. Throughout this paper, we adopted the AB magnitude 
system and WMAP cosmology (SIm= 0.27, VLa_= 0.73, Ho= 
71 km s~' Mpc~' ) unless otherwise stated. 

2. DATA 

This work is based on optical U,g,r,i,z, near-infrared (NIR) 
J,H,K, IRAC 3.6^,4.5^ and GALEX near-UV images of a 
patch 0.6 X 0.6 square degrees wide in the SWIRE Lock- 
man Hole field, hereafter the Deep Swire Field (DSF). This 
patch is approximately centered on the region covered by deep 
VLA imaging ( 10''46'"00^ 59°01'00"). An extensive spec- 
troscopic campaign secured spectroscopic redshifts for sev- 
eral hundreds objects, as detailed in paper III. 

Optical U, g, r images were obtained in 2002 (g,r), and 
2004 (U) at the Kitt Peak National Observatory (KPNO) May- 
all 4m telescope. A detailed description of these images, 
including dat a acquisition and processing, can be found in 
iPoUetta et al.l (12006) . A deep i-band image was obtained 
from the Canada France Hawaii Telescope (CFHT) Mega- 
Cam Science Archive. The data were acquired in 2005 dur- 
ing the observing run 2005BH99. The stacked MegaCam 



image has been produced by the MegaPipe pipeline at the 
CADC''. Medium deep K band imaging, covering about 90% 
of our field, has been downloaded from the UK Infrared 
Telescope (UK I RT) I nfrared Deep Sky Survey (UKIDSS, 
Lawre nce et al.l (120071) ) science archive . UKIDSS uses the 
UKIRT Wide Field Camera (WFCAM, Casali et al. (2007)) 
and a photometric system described in Hewettet al. (2006). 
The pipeline process ing and science archive are described in 
iHamblv et al.l (120081) and Irwin et al., in preparation. We have 
used data from the DR2 data release, which is described in 
Warren et al. (2007). IRAC 3.6/i and 4.5/i images are part of 
the Spitz er Wide-area In fraRed Extragalactic Legacy Survey 
(SWIRE, i Lonsdale et all G2003 2004)). GALEX NUV deep 
imaging has been acquired on the DSF as part of the Deep 
Imaging Survey (DIS). Two contiguous GALEX NUV point- 
ings overlap on the DSF field. Deep stacks have been made 
publicly available in early 2008 with the GALEX Release 4 
(GR4) " . In order to produce a single image covering the DSF, 
the t wo GR4 imag es have been coadded using the Swarp soft- 
ware (' Bertinll20Cl3i) . 

Finally, we used proprietary and/or still unpublished z, J 
and H data. Imaging in z-band was obtained with the MO- 
SAIC camera on the 4m telescope at KPNO on the nights of 
April 2-5, 2005. Ten hours of on-sky integration was obtained 
for a pattern of pointings which produced an image 48' x 48' 
in size centered on the field. The image was reduced with 
the standard IRAF MOSAIC package. Imaging in J and H 
band was obtained with WFCAM on UKIRT on the nights of 
April 6-9, 2007. A total of eight hours on-sky was obtained 
for each band to construct a mosaic image covering 54' x 54' 
centered on the field. The data were pre-processed with the 
UKIRT summit pipeline and then shipped to Cambridge for 
further processing, including removal of instrumental signa- 
ture, sky subtraction, and stacking of microstepped images'^. 
The Cambridge processed images were then mosaicked with 
the SIMPLE Imaging and Mosaicking Pipeline'-^. 

3. CATALOGS 

Cat alogs were generated with Sextractor (i Bertin & ArnoutsI 
[19961) in "dual -image" mode. Three detection images were 
used, hereafter referred to as "optical" (g+i), "NIR" (Jh-Hh-K), 
and "IRAC" (3.6yLtmH-4.5/xm). Even though the best strategy 
for building a detec tion image w ould be to use all images at 
the same time (e.g.. lSzalav et al.|[T999) . the significant differ- 
ence in quality and resolution among our images suggest that 
we can build a better detection image by using only the best 
quality images we have (in terms of resolution, artifacts, bad 
areas due to bright objects). For this reason, the primary de- 
tection image was built from the g and i images alone. In 
order to include in our analysis also redder objects, we then 
considered the NIR and IRAC detection images as well. 

All optical and NIR images were convolved with a gaussian 
kernel in order to match the seeing of the worse optical image, 
namely the U band with a seeing of ~ 1 .3". In order to mea- 
sure photometry in SExtractor dual-image mode, all images 
were registered on the optical detection image pixels, with 
the IRAF tasks GEOMAP and GEOTRAN. The accuracy of the 
pixel registering, as measured with bright point-like sources, 

' For a detailed description of the MegaPipe processing see 
http://www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/megapipe/index.html 
http://www.ukidss.org/archive/archive.html 
" http://galex.stsci.edu/GR4/ 

see http://casu.ast.cam.ac.uk/surveys-projects/wfcam/technicall 
see http://www.asiaa.sinica.edu.tw/~whwang/idl/SlMPLE 
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goes from 0.2-0.3 pixel typical for the optical and NIR im- 
ages, to « 0.5 pixel for the IRAC images (corresponding to 
less than 0.1 native IRAC pixel). 

Photometry was extracted in circular apertures of 0.75" and 
1" radius. Since optical and NIR images were all smoothed at 
the same resolution, their photometry needs no correction for 
different FWHM of the original images. This is not the case 
for the GALEX and IRAC images, which have a significantly 
worse resolution. In order to avoid greatly degrading the reso- 
lution of the optical/NIR images, GALEX and IRAC photom- 
etry was corrected separately for the effect of the larger PSF, 
by means of aperture corrections based on point-like source 
profiles matching the individual GALEX/IRAC resolution to 
the optical/NIR matched resolution of ~ 1 .3". 

Mini-background maps produced by SExtractor were used 
to automatically identify bad areas in each of the three de- 
tection images, i.e. areas affected by large bright star ha- 
los, spikes and other kind of artifacts, image corners poorly 
exposed, and generally all areas whose background is not 
uniform with the rest of the image. While these areas are 
discarded when dealing with statistical studies of the opti- 
cal/NIR/IR galaxy populations, we still consider them for the 
study of the radio-selected galaxy sample, yet flagging the ob- 
jects falling in these areas for a subsequent visual follow-up. 
Furthermore, all the single-band images were similarly ana- 
lyzed in order to flag areas where photometric quality cannot 
be considered uniform with the rest of the image. Photom- 
etry in these areas was then handled with particular care, as 
explained later. 

The accuracy of the aperture photometry measured on these 
images in dual-image mode was tested with simulations, by 
adding to the images point-like sources and attempting to 
recover them given their position is known. It is impor- 
tant to stress that the image quality/depth measured in this 
way doesn't take into account the detection problem of faint 
sources in the single image, since we use the detections com- 
ing from one of the three detection images described above. 
Artificial sources with a gaussian profile of 1.35" for the op- 
tical and NIR images, and 1.8" for the IRAC images were 
added with the IRAF task MKOBJECTS. Assuming their po- 
sitions were known, their photometry was extracted in dual- 
image mode as for the real objects, by creating bright enough 
objects at the same position in the detection image. Again we 
note that this test only aims at checking photometric accuracy 
of detected sources, independent of detection issues: detec- 
tion is assumed to be made on a higher S/N detection image, 
possibly in a very different passband. Nonetheless, these sim- 
ulations include the effect of bad areas in both the detection 
and photometry images, which affect the detection (in the de- 
tection image) or the flux measurement (in the photometry im- 
age). The retrieved aperture magnitudes were then compared 
with the input (total) magnitudes (by applying the appropriate 
aperture corrections). From these simulations, we estimated 
for each image the following quantities, as a function of in- 
put magnitude: i) the percentage of input objects for which 
SExtractor was able to measure a magnitude, ii) the median 
difference Amag between input and output magnitudes, and 
iii) the 16'''-84''' percentiles of the Amag distribution. 

We used these quantities to determine, for each image, the 
magnitude range where the measured flux can be considered 
meaningful, and a realistic error on such flux. In particular, we 
discarded all measurements fainter than the magnitude magf,,, 
defined as the faintest magnitude where more than 90% of the 
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Fig. 1. — ARA vs ADec for radio-optical/NIR/IRAC matching IDs 
(upper-left and lower panels). All matchings are shown for each catalog 
within 2". In each panel the red points show the median and the errorbars 
show the 18''784''' percentiles of the ARA and ADec distributions. The 
upper-right panel shows the histogram of the distances between radio IDs 
and matched counterparts, independently for each of the three catalogs. 

input objects has a measured flux, the median Amag is less 
than 0.2 mag, and both the 16''' and 84''' percentiles of the 
Amag distribution are within one magnitude from the input 
magnitude. While these criteria allow to retain the advantages 
of the dual-image extraction in measuring faint fluxes, they 
also allow to define for each passband a limit beyond which 
such fluxes are no longer deemed meaningful measurements, 
and thus are treated as drop-outs. From the error curve de- 
rived from these simulations we also determine a magnitude 
mag 10 where simulated objects have their flux measured with 
a typical error of less than 0.1 mag. The adopted values of 
magcu, and magio are given in Table [T] 

4. THE RADIO-SELECTED SAMPLE 

The sample of radio sources used here is described in detail 
in paper I. The DSF was observed with the VLA in the A, 
B, C, and D configurations for a total of almost 140 hours 
on-source. The final image has a resolution of 1.6", with a 
typical rms at the center of the image of 2.7 /iJy ( see paper I 
for details). 

The sky region included in our detection images contains 
1930 sources out of the 2055 of the original paper I radio cata- 
log. All three (optical/NIR/IRAC selected) catalogs described 
in section[3]were used to identify a (optical/NIR/IR) counter- 
part for the radio sources in the catalog. The radio and opti- 
cal/NIR/IRAC WCS were first matched by correcting for the 
mean offsets in ARAop„c„/_rfld/o and ADeCop„cfl/-ra</w- The 
scatter in l\RAopticcd~mdio and ADeCop„-,.a/_™/io is ~0.3"(see 
figure[T]l. Then, for each radio source the three catalogs were 
searched for a match within 1" (with an order of preference 
1) optical, 2) NIR, and 3) IRAC. Most of the radio sample 
was matched with an optical counterpart: out of the whole 
sample of matched counterparts, the (optical unmatched) NIR 
counterparts and (optical and NIR unmatched) IRAC coun- 
terparts contribute for about 12% and 9%, respectively. How- 
ever, about 20% of these IRAC counterparts and 40% of these 
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NIR counterparts are located in "bad areas" of the optical 
detection image, where detection was hampered by bright 
sources or artifacts. Also, when restricting to the redshift 
range 0.3 < z < 1.3 which will be the main focus in the fol- 
lowing, IRAC selected counterparts contribute for just 1%, 
and NIR selected counterparts for 7% (half of which in bad 
areas of the optical detection image). 

While for most of the sources a counterpart closer than 
1" was found, for almost 6% of the 1930 sources it was not 
possible to identify a counterpart within 1". Half of these 
could be matched with a counterpart increasing the match- 
ing radius up to 1.5". The other half of the sources without 
counterpart were automatically and visually inspected: some 
of these sources are indeed found in image areas affected by 
bright objects halos/artifacts, some others are blended with 
another bright object - clearly in such cases the possible coun- 
terpart may likely go undetected. However, for 13 (out of 59) 
unmatched sources it seems that there is actually no counter- 
part in our images (and no apparent issues in the images which 
could hamper its detection). 

In conclusion, 94% of the radio sources contained in our 
detection region were assigned a reliable counterpart (within 
1"), while 3% were assigned a counterpart between 1" and 
1.5" and 3% were not assigned any counterpart. 

Since the resolution of the radio image is similar to that of 
the optical /NIR images, w e did not apply the likelihood ratio 
technique (lRichteiill975h which is commonly used to evalu- 
ate the reliability of each identification based on the source- 
counterpart distance and counterpart magnitude. However, 
we used the distance and magnitude distribution of our identi- 
fied counterparts to evaluate the contamination of our catalog 
from false counterparts, estimating the probability of false as- 
sociation by randomly shifting the radio source coordinates 
and repeating the association process. Given the characteris- 
tics of most of our associations (75% (97%) are at a distance 
of less than 0.5"(1") from the source, 70% brighter than i=25 
within an aperture of 1.5") the overall contamination of the 
whole catalog is estimated to be lower than 2%. The contam- 
ination from false counterparts of the 0.3 < z < 1.3 sample 
which will be used in most part of this work is estimated to be 
negUgible (< 1%). 

Out of the sources with an assigned counterpart, 76% were 
detected (i.e. have a measured aperture magnitude brighter 
than magcr) in the U band (and up to 80% including objects 
possibly undetected because in flagged areas of the U band 
image), 83% (up to 89%) in the g band, 80% (89%) in the r 
band, 74% (91%) in the i band, 75% (92%) in the z band, 85% 
(91%) in the J band, 90% (95%) in the H band, 81% (87%) 
in the K band, 88% (100%) in the IRAC 3.6^ band, and 93% 
(100%) in the IRAC 4.5/^ band. We note that we give these 
numbers as an indication of the typical spectral energy dis- 
tributions (SEDs) of the radio-selected sample: the specific 
numbers depend on the different depth and overall quality in 
terms of bad areas of the different images. 

5. PHOTOMETRIC REDSHIFTS 

5.1. Determination of photometric redshifts from 
multi-wavelength NUV to mid-IR photometry 

The optical, NIR and IRAC selected multi-wavelength cat- 
alogs described above were used to estimate photometric red- 
shifts (photo-zs) by means of comparison with a library of 
galaxy SED templates covering a range of star-formation his- 
tories, ages and dust content. A set of 33 templates were 



used, spanning from a classical local elliptical to several star 
forming galaxies to a QSO dominated template, and all cover- 
ing a rest-frame spectral window [1000-70000], thus ensuring 
an adequate cross-correlation with the available photometric 
coverage. Beside local galaxy temp lates (e.g..|Coleman et aP 
1980, Mannucci et al. 2001, Kin ney et alJ 119961) . a seTof 
semi-empirica l templates based on observations plus fi tted 
model SEDs (iMarastonI [T9981 |Bruzual & Chariot* '2003^ of 
« 300 galaxies in die FORS Deep Field (Heidt et al. 2003]) 
and Hubble Deep Field ( Williamset al.lil99&) . were included 
to better represent objects to higher redshifts. 

Here we describe bri efly the method used to estimate 
photo- zs, a nd we refer tolBender et all (1200 lb . iGabasch et al.l 
d2004b and iBrimioulle et all ( 12008^ for a more detailed de- 
scription of this method, as well as of the construction of the 
templates. The aperture PSF-matched photometry of each ob- 
ject in the available (up to 1 1) passbands was compared with 
the templates, calculating a redshift probability function over 
the range < z < 10 (in steps of 0.02) for afl SEDs. This is 
done assuming some priors: a different prior on the redshift 
distribution is assumed for different types of templates (cor- 
responding to younger or older SEDs, that is e.g. an old lo- 
cal elliptical template is assumed to be increasingly unlikely 
at higher redshifts, while templates corresponding to young 
stellar populations or QSOs are assumed to have a basically 
flat likelihood across all redshifts explored). Furthermore, a 
weak, broad prior on the absolute optical and NIR magnitude 
lowers the probability to have magnitudes brighter than -25 
and fainter than -13. Ly-a forest depletion of galaxy templates 
is implemented according to Madau (1995). The "best-fit" 
photo-z Zpitoi is chosen as the redshift maximizing the prob- 
ability among all templates, and an error on Zphoi is defined 

as e^phot = [^ijizi — ZphotyPijV^^, with zi the considered red- 
shift steps in < z < 10, and Pjj the contribution of the j-th 
template to the total probability function at redshift z,. 

Systematic offsets between the measured and predicted col- 
ors as a function of redshift, which may be due to several rea- 
sons as for instance uncertainties in the estimated zero-point, 
but also possibly uncertainties in the filter response curves or 
even systematics in the templates, were estimated using al- 
most 500 spectroscopic redshifts available in our field (paper 
III, private communication from G. Smith 2007, plus some 
few more published redshifts available from NED'"*). There- 
fore, the zero-point for each of the passbands was corrected 
by a factor which minimizes the systematic shift between ob- 
served and predicted color for well-fitted spectroscopic galax- 
ies. 

Also stellar templates (iPicklesI 119981) were fitted to all 
sources. Comparing stellar and galaxy for objects with a 
point-like morphology, we checked that objects brighter than 
I~23.7 with a best-fit stellar lower than the best-fit galaxy 
could be classified as stars. The number counts of such 
selected "stars" are in good agreement with predictions of the 
iRobin et all (12003 1) Galaxy model Fainter than I~23.7, the 
classification was found less reliable. None of the sources in 
the radio sample was classified as star. 

When fitting an object's photometry, besides excluding 
from the fit all magnitudes deemed unreliable (including mag- 
nitudes of objects in flagged areas), all magnitudes fainter 
than magcu, (as defined above) were considered as drop-outs, 

'** http://nedwww.ipac.caltech.edu/ 
15 www.obs-besancon.fr/model/ 
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Fig. 2. — The comparison of spectroscopic and photometric redshifts for 
more than 400 sources in the field. In the upper-left and lower panels the 
estimated photometric redshift is plotted against the spectroscopic redshift, 
sepai'ately for the three optical, NIR and IRAC selected catalogs (the numbers 
of matched sources are 402, 419 and 429 respectively). The solid line in each 
panel marks the bisector. Orange symbols highlight radio sources, making up 
about 3/4 of the spectroscopic sample. The histograms of Az/(1 -f z) for all 
three catalogs are plotted in the upper-right panel. 

and a few photo-z determinations with different ways of deal- 
ing with drop-outs were compared, in order to test the stabil- 
ity of the derived photo-z. Obviously, the estimated photo- 
z is most unstable for those objects with a very high num- 
ber (e.g. > 6) of drop outs (or combination of drop outs and 
flagged magnitudes). Only a minor fraction of the objects is 
concerned, and the effect on the global sample, and in par- 
ticular at redshift lower than 2, can be considered negligible. 
Nonetheless, the outcome of this test was used together with 
the of the best fit, the error on the estimated photo-z, and 
on the total number of actually measured magnitudes used 
(e.g. not upper limits or flagged magnitudes), to evaluate the 
quality of the estimated photo-z for each source. 

Only "constrained photo-zs", i.e. with at least 4 measured 
magnitudes actually used (thus not including upper limits and 
flagged magnitudes) and a discrepancy between different de- 
terminations of less than 20% in Az/{l+z) were deemed re- 
liable'^ In the following we only use photo-zs deemed 
reliable based on these criteria, unless otherwise noted. 



A quality flag QF was defined as follows: QF=AA: a photo-z is de- 
termined in all different realizations, the total number of upper limits plus 
flagged magnitudes is less than 5, the maximum difference Az/(l+z) among 
the different determinations is less than 20%, the estimated error on the 
photo-z is less than 0.2(l+z). QF=A: same as QF=AA but the constrain 
on the estimated eiTor is relaxed to 0.4(l+z). QF=B: same as QF=A but con- 
strain on number of flagged magnitudes plus upper limits is increased to less 
than 7. QF=C1: same as QF=B but no constrain on the estimated error on 
photo-z. 

Also, based on the distribution of such selected photo-zs and on the 
comparison with the spectroscopic sample (see below), all photo-zs having 
at the same time both a > 2 (97''' percentile of the distribution of the 
"constrained photo-zs") and a relative error e-pi,„, /zphot larger than 75% were 
discarded. 



5.2. Accuracy of photometric redshifts 

The overall accuracy of our photo-zs is estimated by com- 
parison with the spectroscopic sample. Figure |2] shows the 
comparison of the spectroscopic redshifts with the retained 
reliable photo-zs for each of the photometric catalogs (op- 
tical, NIR and IRAC selected). In each case, more than 
400 objects have been used for the comparison, deriving 
a median(Az/(l -I- z))< 0.003, «3.5% outliers (defined as 
Az/(1 +z)> 0.2), and an accuracy of w 0.05 in Az/(1 +z), 
as estimated either from the NMAD estimator (iHoaglin et al.l 
ll983LIIlbertet all 120091) or from the 16"'-84"' percentiles of 
Az/(l+z). 

While we have a sizable spectroscopic sample allowing 
us to assess the accuracy of our photo-z determination, it 
may also be useful to compare our results with independently 
determined photo-zs. While such comparison has not the 
same strengths of the usual spectroscopic vs photometric red- 
shift comparison, it may help to overcome two of its main 
weaknesses: spectroscopic samples are typically only a small 
fraction of a galaxy sample, and are significantly biased to- 
ward bright sources. In figure |3] we show the correlation be- 
tween the photo-zs der i ved in this work and those derived by 
iRowan-Robinson et all (120081) (hereafter RR08). RR08 de- 
rived photo-zs for all the SWIRE survey, including the Lock- 
man Hole field used in this work. They not only use a different 
code for estimating photo-zs , but also a significantly different 
approach. While having a full optical-i-NIR coverage in some 
of the fields, in the Lockman Hole they only use photometry in 
U, g, r, i, 3.6 and 4.5/im passbands. They report for the whole 
SWIRE survey a typical rms of {Zphor'^si>ec)l{^+'^siyec), exclud- 
ing outliers, slightly larger than 4% for 6 passbands (which is 
the case of the Lockman Hole) and rvega <24, and w 4% out- 
liers. In the left panel of figure |3] we thus highlight the com- 
parison for sources with Xyega <24 and photometry available in 
6 passbands in the RR08 catalog. By comparing our photo-zs 
with those by RR08 for these objects, we find almost 6% out- 
liers, a median (zRR08-z,/,/,vvo,vt)/(lH-z,ft,-.„„o,.i) of « -0.03, and a 
scatter Ri6%, which is consistent with what expected based on 
the claimed scatter of both works. We note that their perfor- 
mance, for all SWIRE fields, worsens for z > 1 .5 (see orig- 
inal paper). In fact, this worsening may be expected to be 
especially significant for the Lockman Hole where no NIR 
photometry is used in RR08, thus likely hampering photo- 
zs at z > 1 — 1.5. In fact, as figure |3] shows, the compari- 
son between our and RR08 photo-zs gets significantly worse 
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Fig. 3. — The comparison of this work photo-zs with those derived by 
RR08. Left panel: gray symbols show all data points while black symbols 
show those which in the RR08 catalog have photometry in 6 passbands and 
mag r< 24. Right panel: gray symbols show all data points while black 
symbols show objects with available spectroscopic redshift in either this work 
or RR08. See text for details. 
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Fig. 4. — The compai'ison of this work photo-zs with spectroscopic red- 
shifts of 198 X-ray sources from the Trouille et al. (2008) sample. Solid 
symbols refer to a subsample of 7 1 radio detections, while empty symbols do 
not belong to the radio sample. The solid line is the bisector. Color coding 
reflects the Tr ouille et al., t2008i) optical spectral classification as indicated 
(see text for details). 

beyond zw 1.5. As far as our photo-zs are concerned, for 
« 15 objects with 1.5 <Zspec < 3.2 we have a median {Zphor 
'ZspecVi^+'^spec) ^ 0.01 and scatter < 0.05, pretty similar to the 
figures for the whole sample, but a formally higher number of 
outliers (<7%, due to one outlier out of ?sl5 objects). The 
right panel of figure |3] highlights instead the comparison of 
RR08 and our photo-zs for sources with spectroscopic red- 
shift available, either in our catalogs or in RR08. For these 
sources, while finding a scatter similar to that of the whole 
magnitude- and number of passbands- limited sample, we 
find a lower median (w —0.01), and < 5% outliers. We no- 
tice that the spectroscopic samples used in this work and in 
RR08 have a significant overlap, thus this might affect the 
comparison of the photometric redshifts results. Nonetheless, 
the comparison of photometric redshifts derived with differ- 
ent photometry, different SED libraries, and different meth- 
ods, allows an indirect evaluation of the photo-z performance 
on a much larger sample, and at fainter magnitudes, than those 
allowed by the usual comparison with spectroscopic redshifts. 
This comparison confirms, at least out to redshift ~ 1.5, the 
overall accuracy of our photo-zs as estimated by comparison 
with our spectroscopic sample. 

On the other hand, it is interesting to compare our photo- 
zs with a spectroscopic sample of X-ray sources in this field 
published by Trouille et al. (2008). In figure |4] we compare 
our photo-zs with (Trouille et aLi .2008 ) spectroscopic red- 
shifts for «200 X-ray sources (out of which w 70 belong 
to our radio sample) for which we have a deemed reliable 
photo-z''^. We divide and color-code the sample according 

Five more objects common to both the I Trouille et al.l 120081) and our 
sample were excluded from this comparison because of possibly dubious 
spectroscopic redshift. One of these is a z=0.35 object for which an available 
spectroscopic redshift is consistent with our photo-z. The others all have 
redshifts z> 3 in the Trouille et al. catalog, and thus might be expected to be 



to the ITrouille et al.l (l2008h optical spectral classification, in 
absorbers, star-formers, high-excitation sources and broad- 
line AGNs (see the original paper for details). The frac- 
tion of such "peculiar" objects in this comparison sample is 
quite relevant, with 40% of the sample being made of broad- 
line AGNs, and a further « 30% of high-excitation sources. 
As figure |4] shows, our photo-zs perform significantly worse 
for this spectroscopic sample, as compared to our (or RR08) 
spectroscopic sample, with much larger scatter and number 
of outliers. The overall statistics for Az/(1h-z) for this sam- 
ple is median< -0.015, NMAD scatter «10%, and «20% 
of the objects having Az/(1h-z)>20%. As it is clear from 
the figure, the worse results are obtained for the broad-line 
AGN sub-sample (NMAD scatter «17%, ?s30% of the sam- 
ple with Az/(1h-z)>20%). The poor agreement obtained from 
this comparison is not unexpected, due to the contamination 
of the broad-band photometry wit h a strong AG N contri- 
bution (see e.g. Polletta et a l. 2007, ISalvato et al.' 2009, and 
references therein). We also note that our photo-zs perfor- 
mance appears to be similar or better than that achieved by 
ITrouille et al.] (l2008h . even though when estimating photo- 
zs they include templates built from their spectroscopically 
confirmed broad-line AGNs, while our template set was not 
specifically tailored toward AGN-dominated SEDs. Further- 
more, we note that many (not all though) of the outliers have 
a large error associated with the estimated photo-z. 

5.3. Photometric redshifts for the radio sample 

An estimated reliable photo-z is determined for 1610 radio 
sources (86% of the identified counterparts and 83% of the 
whole radio sample included in the detection image), while 
for the rest either no photo-z could be estimated, or in most 
cases it was deemed not reliable according to the criteria de- 
fined above. 

In figure |2] the s» 300 radio sources making up about 3/4 
of the spectroscopic sample are highlighted with orange sym- 
bols. While one might expect a worse performance of the 
photo-zs for the radio sample, because of the presence of the 
significant AGN population among radio sources, the statis- 
tics for the radio sample alone are quite similar, and only 
slightly worse, than those obtained using the whole spectro- 
scopic sample: median Az/(1 +z)=0.0008, scatter Ri5.5%, 
and 4% outliers (however we note again that a large fraction 
of the whole spectroscopic sample is made of radio sources). 
In figure |4] we have already shown (solid symbols) how our 
radio sample photo-zs perform for the Trouille et al. (2008]) 
X-ray selected sources. In the following, we include the spec- 
troscopic redshifts from Trouille et al. (2008) in our analysis 
of the radio sample. 

In figure|5]we show the redshift completeness of our sample 
(including spectroscopic and reliable photo-zs) as a function 
of redshift'''. Our criteria for selecting a reliable photo-z nat- 
urally disfavor higher redshift sources. Since of course we do 
not know the redshift of all the sources, we can do just an ap- 
proximate estimate: we assumed that all our photo-zs were 
broadly correct, even those we decided were unreliable, and 
used them as a reference to estimate the redshift complete- 
ness of the sample. By comparison with our different photo-z 
catalogs based on different settings (of which figure |5] is one 

very faint or drop-out in the U band. However they are all clearly detected in 
our U band image, and in some cases in the GALEX NUV as well. 

" In the following, 8 objects with QF=C1, unusually low x' and very high 
estimated photo-z error were excluded from the reliable photo-z sample. 
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Fig. 5. — The redshift completeness of the radio sample. Upper panel: 
The shaded histogram shows one realization of the redshift distribution of 
the whole sample of identified counterparts, while the solid line shows the 
histogram of used redshifts (spectroscopic redshifts or reliable photo-zs). 
Lower panel: The fraction of available redshifts (spectroscopic redshifts or 
reliable photo-zs) as a function of redshift: the sample is assumed to be more 
than 90% complete in the redshift range 0.3 < z < 1.3. 



example), we estimated that in the redshift range 0.3-1.3 our 
redshift completeness should always be higher than 90%, and 
the 0.3 < z < 1.3 sample is overall complete at a 95% level. 
In other words, if our photo-zs are broadly correct so that the 
total number of objects in the redshift range 0.3 < z < 1.3 
is right, we estimate an overall redshift completeness in the 
0.3 < z < 1.3 redshift range of more than 95% of the iden- 
tified counterparts, and more than 90% of the whole radio 
sample within the detection image^" assuming that all radio 
sources without an identified counterpart are in this redshift 
range. This is definitely a conservative assumption since it is 
likely that many of the unidentified counterparts are very faint 
galaxies probably at higher redshift. 

5.4. Redshift distribution of the radio sample 

In figure|6]we show the redshift distribution of radio sources 
based on the radio (sub-)sample for which either a reliable 
photo-z or a spectroscopic redshift is available (gray shaded 
histogram). The redshift distribution is shown up to redshift 
z ^ 2.7, beyond which the redshift completeness is expected 
to drop below 80%^'. The redshift distribution plotted is not 
corrected for the estimated redshift incompleteness as esti- 
mated above (e.g. figure|5]|. 

Within the detection image here means within the boundaries of the 
detection images and excluding regions within the boundaries where all three 
detection images were masked. 

^ ' We also remind the reader that the sample shown here includes a negli- 
gible fraction (less than 2%) of counterparts matched within 1.5". 



Since the DSF radio image is obtained with a single VLA 
pointing, the rms is not constant over the field, but increases 
with the distance from the field center. For this reason, the 
whole radio sample does not have a single flux density limit 
and this needs to be taken into account. In the following, 
wherever it is needed we will correct the biased nature of the 
whole sample by taking into account the non-uniform rms. 
While it is virtually impossible to make a perfect correction, 
due to the unknown size distribution of the radio sources, for 
our purpose we will calculate the appropriate correction for 
all resolved sources for which a reliable size could be esti- 
mated, and will assume that all remaining sources have a typ- 
ical source size of 1.2"(see paper I). Taking into account the 
bandwidth and time smearing, as well as the primary beam 
correction, and the resolution of the image where the detection 
was performed, we can estimate how the rms changes across 
the image, and thus the 5-a limiting flux density^^ at each dis- 
tance from the field center (see paper I for more details). In 
figure |2] we show, as an example, the estimated 5-a limiting 
flux density as a function of distance from the field center for 
a source of size~1.2". Based on this estimated 5-0" limiting 
flux density as a function of distance, and on the masked areas 
in the optical/IR detection images, we calculated the area A,- 
over which each object could be observed and enter our sam- 
ple. The sky coverage is then used to weight each object (by 
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Fig. 6. — Top panel: the gray shaded histogram shows the redshift dis- 
tribution in the range 0.1 < z < 2.7 based on the sample of radio sources 
with identified counterpart and reliable redshift estimate (see text). The red 
(dotted) and purple (continuous) lines show the redshift distribution of flux- 
limited samples with radio flux density above 24/^Jy and 84/^Jy, respectively. 
All distributions are sky-coverage corrected (see text). All histograms are 
normalized to the total number of objects in each sample. Bottom panel: 
The red (continuous), yellow (dotted) and blue (dashed) lines show the red- 
shift distributions of sub-samples with radio flux densities in the range [16- 
30], [30-60], and > 100/iJy, respectively. In both panels, all distributions are 
sky-coverage corrected (see text), and are normalized to the total number of 
objects in each sample. Errorbars show Poissonian errors. 

^- The radio sample used in this work only includes detections with S/N > 
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Fig. 7. — The solid line shows the 5-cr limiting flux density as a function 
of distance from the center of the VLA field, for a source size of 1.2" (see 
text for details). Empty symbols show the integrated flux density vs. distance 
from the field center for actual detections. Thick dotted lines of different 
colors show four examples of flux-limited sub-samples selected in uniformly 
covered portions of the image (R < 5' and Siaghz > Ifi/^Jy- ^ < 1"' ^^nd 
Si 4G«- > 24/iJy, R < 15' and Si 4g//' > 41/iJy, and R < 20' and Si 4GHz > 
84^(Jy). 

IM,) and thus estimate properties for a flux-limited sample. 
We note that, as it can easily be seen from figure|2l the faintest 
sources in our sample can be observed over a very small frac- 
tion of the survey (<10%). One might thus be concerned 
that inaccuracies in the sky coverage correction would sig- 
nificantly affect the results. However, we note that all results 
presented in the following would be unchanged if considering 
only flux-limited sub-samples selected in uniformly covered 
portions of the image (e.g., R < 5' and Siaghz > 16/iJy, or 
R < 20' and Siachz > 84^Jy, see figl?]), while the sky cover- 
age correction allows us to make full use of our data set. In 
the following, we will explicitly refer to flux-limited samples 
when using sky-coverage corrected samples, in contrast with 
the whole, biased sample. 

All distributions plotted in figure |6] are corrected by the 
sky coverage. Beside the redshift distribution derived from 
the whole sample, in the top panel of figure |6] we also show 
the redshift distributions of two flux-limited samples with 
S\ACHz >24 and >84 /iJy. These suggest that the redshift 
distribution might depend of the flux density of the sample, 
and we try to make this clearer in the lower panel of figure|6l 
where we plot the redshift distributions of three sub-samples 
in different flux density ranges, our faintest sources (16. 5< 
SiACHz <30 /iJy), our typical about-median flux population 
(30< SiACHz <60 /iJy), and bright sources (Siachz >100 
/iJy). The redshift distribution of the faintest sources seems 
to be different from that of the bright ones, with a sharper 
peak at z « 0.9. A Kolmogorov-Smirnov test indeed suggests 
that the two distributions are different at a high significance 
level (P < 0.005)23. 

6. SPECTRAL ENERGY DISTRIBUTION FITTING 

We use SED fitting on the available multi-wavelength pho- 
tometry in order to estimate fundamental properties of the 
stellar populations hosted in this radio source sample. In the 

We note that the Kolmogorov-Smirnov results quoted here and in the 
following are obtained from flux-limited sub-samples selected in uniformly 
covered portions of the image, since the test cannot be meaningfully applied 
to the sky-coverage corrected data 



following we will consider photo-zs just as spectroscopic red- 
shifts, assuming that the galaxy is at that redshift and ignor- 
ing any error on the photo-z. Different SED fits are used in 
the following. The first characterization of each galaxy SED 
is given by the best-fit template associated with the best-fit 
photo-z. As already said above, and as will be described more 
in detail below, the templates span a range in stellar popula- 
tion ages, including dust extinction. These templates will then 
be used to classify galaxies according to the broad, average 
properties of their stellar populations. In order to properly 
treat galaxies with an available spectroscopic redshift, the fit 
was recomputed for these objects assuming the spectroscopic 
redshift. 

In order to classify SEDs based on a more "parametric" 
approach, and also to investigate possible misinterpretations 
coming from the use of non-evolving templates, we also per- 
formed SED fits using stellar population synthesis models 
produced with the Br uzual & Chariot (2003) code. Star for- 
mation histories (SFHs) were parametrized by simple expo- 
nentially declining star formation rates, with a timescale r 
ranging between 0.1 and 20 Gyr, and age between 0.01 Gyr 
and the age of the Universe at the object's redshift. Metallic- 
ity is fixed to solar and a Salpeter ( 1955) initial mass function 
is adopted, with lower and upper mass cutoffs of 0.1 and 100 
Mq. a variable amount of extinction by dust is also included, 
with Ay g [0, 1.51. The fitting pr ocedure is described in full 
detail in lDrorvet al] (I2004ll2005h . These results too wiU be 
used in the following to roughly characterize the host stellar 
populations based on the best-fit age/r (meaning the age of 
the stellar populations from the onset of the star formation 
divided by the e-folding time of the exponentially declining 
SFH). 

Finally, SED fitting is also used to estimate stellar masses. 
For this purpose we use, as it is customary, a two-component 
model, adding to the main smooth component (exponentially 
declining SFHs described above) a secondary burst. This 
burst is modeled as a 100 Myr old constant star formation 
rate episode. Metallicity and IMF of the burst are the same 
assumed for the main component, however dust extinction for 
the burst is allowed to reach higher values (Ay G [0,2]). We 
notice that since we fit aperture photometry, the stellar masses 
obtained refer to the portion of the galaxy contained in the 
aperture. We correct these masses to "total stellar masses" 
by means of the ratio of "total" (FLUX_AUTO) and aperture 
fluxes in the detection images. The median correction applied 
for the radio sample is a factor w 2, and 90% of the sample 
is corrected by a factor ranging between 1.5 and 4. This is an 
approximation neglecting any color gradients which may cer- 
tainly exist in the galaxy, in other words we assume that the 
stellar populations within the measured aperture are represen- 
tative of the stellar populations of the whole galaxy. 

Some of the main derived properties which are used in this 
work, for objects with spectroscopic or reliable photo-z, are 
listed in Table |2] which is available in full on the on-line ver- 
sion. 

7. SED PROPERTIES OF HOST GALAXIES: AGN ACTIVITY AND 
STAR FORMATION 

In the following we compare the radio and optical/NIR 
properties of this radio sample based on the SED fitting re- 
sults described above. 

As shown in figure [8] galaxies best-fitted by different kind 
of templates (thus in principle galaxies with different stel- 
lar populations) tend to occupy different locations in the ra- 
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Fig. 8. — The 1.4 GHz luminosity plotted against redshift for objects with 
redshift between 0.3 and 1.3. Grey points in all panels show all data points 
while colored symbols show objects whose photometry was best-fitted by 
SEDs of different types, as indicated in each panel. In each panel, connected 
small squares show the (sky coverage corrected) interquartile range of the 
^lAGHz distribution as a function of redshift for the relevant SED-selected 
sub-sample. 

dio luminosity against redshift diagram. The simplest, and 
expected, explanation is that different radio luminosities are 
associated to different physical processes, namely star for- 
mation and AGN activity. Figure |8] shows the radio lu- 
minosity against redshift for classes selected based on the 
best-fit photo-z SED template. The templates were divided 
in "quiescent" ( including for instance the elliptical template 
by IColeman et al. (1980) and the SO and Sa templates by 
[Mannucci et al. 2001), "interme diate" galaxies with low star 
formation (including e.g. the Sb ' Mannucci et al.] ( 1200 Ih tem- 
plate), and "star-forming" templates including all actively 
star-forming galaxies, plus a "QSO" class of a few objects 
best-fitted by a QSO template. As a reference, the restframe 
U-B color^'* ranges for the three classes of quiescent, inter- 
mediate, and star-forming objects are approximately 1.1-1.4, 
0.9-1.1, 0.1-0.9, respectively. Similarly, the U-V color ranges 
are approximately 1.9-2.2, 1.5-1.8, 0.1-1.3, while the break 
strengths at 4000, D„(4000) ^\ are about 1.6-2.1, 1.5-1.7, 1- 
1.3. For reference, these D„(4000) ranges may be compared 
for instance with typical D„(4000) values for different kinds 
of stellar population s measured in a local (SDSS) sample 
dGallazzi et al.l20()5l their figure 7) , and in a VVDS sam ple in 
the redshift range 0.45 < z < 1.2 dFranzetti et al.ll2007l their 
figure 4). We also note that in figure [8] the "star-forming" 
class is further split in two sub-classes of star-forming tem- 
plates, "redder SF" and "bluer SF", with restframe U-B rang- 
ing in 0.7-0.9 and 0.1-0.7, and A, (4000) about «1.3, 1-1.2, 
respectively. 

Due to the non-evolving nature of the templates used in the 
photo-z determination, figure |8] only shows the stellar popu- 

Restframe U-B colors, here and in the following, are calculated in the 
[Buser & Kurucz ( 1978) U and B3 filters. 

We adopt the Balogh et al. 1 1999) definition of the D„(4000) index, that 
is the ratio of the average flux densities in the narrow bands 4000^100 and 
3850-3950 . 



lations status (i.e. actively star-forming, passively evolving, 
etc.) at the time of observations, without any evolutionary 
link between same-class objects at different cosmic epochs. 
In other words, depending on the specific star formation his- 
tory of each galaxy, and on the overall evolution of galaxy 
stellar populations, galaxies may (and will) change their class 
as time goes by. 

As figure [8] shows, at all redshifts the highest radio lumi- 
nosities in the probed range (e.g. I^iaghz 

> 2 X 10^^ W/Hz 

at z 0.5 or I^iaghz > 5 x 10^^ W/Hz at z 1) are typ- 
ical of low-starforming systems (i.e., galaxies classified as 
intermediate or quiescent). Nonetheless, radio luminosities 
of such low-starforming galaxies span all the range covered 
by this survey. On the other hand, galaxies classified as ac- 
tively star-forming tend to avoid the highest radio luminosi- 
ties, with just few exceptions. The median radio luminosi- 
ties of all sub-samples are similar, with star forming galax- 
ies typically showing a median radio luminosity about 20% 
lower than low-starforming systems (quiescent and interme- 
diate), over the redshift range probed. However, the mean 
radio power of intermediate and star forming galaxies is typ- 
ically lower than that of quiescent galaxies (the mean h\ncHz 
of bluer star-forming, redder star-forming, and intermedi- 
ate galaxies are 30%, 40%, and 50% of the mean \^\aghz of 
quiescent classified sources, respectively). The interquartile 
ranges for Li.4gh; plotted in figure[8]also show how the typ- 
ical range and spread in radio power is different for different 
SED classes, with the bluer star-forming galaxies mostly ly- 
ing just above the flux density limit of this survey. 

Finally, figure [8] shows how a significant part of this radio 
sample is made of objects classified as intermediate. This is 
likely due to the depth of this survey which allows us to go 
beyond the AGN dominated population, but at the same time 
does not allow us to reach the still fainter radio luminosities 
typical of the bulk of normal star-forming galaxies. We note 
that, while the statistics given above take into account the sky 
coverage, the points plotted in figure [8] show the whole radio 
sample with 0.3 < z < 1.3, regardless of the different sensitiv- 
ity at different radii in the VLA image, therefore the different 
densities of data points are not directly indicative of the ac- 
tual relative fractions of different kinds of objects, which are 
better addressed below. 

7.1. The nature of faint radio sources: a significant 
intermediate population? 

In figure |9l we show the restframe U-B vs B color- 
magnitude diagram of the radio sample, divided as above in 
quiescent, intermediate and star-forming. In this figure, all 
sources with redshift 0.3 < z < 1.3 are plotted^^, taking ad- 
vantage of the small evolution of the restframe U-B color in 
the redshift range probed, and thus just one separation be- 
tween red and blue galaxies was adopted at all redshifts (see 
e.g. Willmer et al. 2006, Cooper et al. 2007, in the same red- 
shift range studied here). The figure shows the well known 
different locations in the color-magnitude diagram of differ- 
ent galaxy types, with quiescent galaxies on the red-sequence 
and star-forming galaxies in the blue-cloud. However, as 
compared to a typical color-magnitude diagram for an opti- 
cally selected galaxy sample, it is evident that the so-called 

However, 22 spectroscopic objects are not plotted because, due to a high 
number of flagged magnitudes, they lack sufficient photometric information 
to perform a reliable SED fit, even with known spectroscopic redshift. This 
negligibly lowers the completeness of the sample plotted, which is now 94%. 
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Fig. 9. — Top panel: The restframe U-B vs B color-magnitude diagram 
for objects with redshift 0.3 < z < 1.3. All objects are included regardless 
of their radio flux density and position in the VLA image. Symbols are color 
coded according to best-fit SED type, as indicated (same broad classes as 
in figure [8). The dotted Une shows the division between red sequence and 
blue cloud from DEEP2 data as in Willmer et al. 1 2006), Cooper et al. 1 2007). 
Bottom panels: The distribution of restframe color U-B for the different SED 
classes, in four flux-limited sub-samples, as indicated (see text for details). 
The dotted lines mark the U-B color of the division line at the median B 
magnitude of each plotted sample. In each panel, histograms for the different 
SED classes are color-coded as in the upper figure (quiescent = horizontal, 
intermediate = diagonal and star-forming = vertical hatched), and the gray- 
shaded histogram shows the distribution of the whole subsample plotted in 
the panel. 



"green valley", often thought to be populated by transitional 
objects which are shutting off their star formation and migrat- 
ing to the red sequence, is not an underpopulated region in 
this diagram, and actually contains a significant fraction of 
this radio-selected sample. 

However, the actual nature of this significant popula- 
tion of apparently intermediate-type galaxies with low star- 
formation activity is not necessarily clear In fact, while it 
might be tempting to consider the option that these are indeed 
"composite" or "transition" systems, meaning galaxies shut- 
ting off their star formation because of (or linked to) an host 
AGN, on the other hand the colors characteristic of an "inter- 
mediate" SED can also be produced in other ways, the most 
obvious being a starburst affected by a very high dust extinc- 
tion which is not properly handled by our template set, or pos- 
sibly a red-sequence galaxy whose observed SED is severely 
contaminated by AGN emission. In fact, in an attempt to look 
a bit further into the dust attenuation issue of the intermedi- 



ate galaxies, we can compare our template-based classifica- 
tion with the parametric one based on the single-component 
model SED fitting. As shown in figure [TOl there is a good 
overall correlation between the age/r of the best-fit model 
and the template-based classification. However, if we look at 
the dust extinction which according to the best-fit is affecting 
the model SED, we find that 10% of the intermediate sample 
at 0.3 < z < 1.3 was best-fitted by a significantly extincted 
young stellar population (age/r < 4, Ay > 1). For compari- 
son, 30% of the galaxies classified as star-forming was best- 
fitted by such kind of model, and 3% of the galaxies classified 
as quiescent. Even though we will not attempt to correct col- 
ors based on the rough attenuation estimated by this simple 
SED fit, the occurrence of apparently significantly extincted 
objects in the intermediate sample may suggest a possibly 
relevant contaminati on by reddened starbursts (see also e.g. 
ICowie & Ba rger 2008). The actual nature of these intermedi- 
ate sources will be further investigated in a forthcoming paper, 
by exploiting X-ray and infrared data. 

For the time being, to avoid misinterpretations, in the fol- 
lowing we will drop the labeling of the three classes as qui- 
escent, intermediate, and star-forming, which directly refers 
to the inferred nature of the host stellar populations, and 
will adopt a more generic and observationally-based "red", 
"green" and "blue". 

It is worth noticing how the nature of the radio-galaxies 
host stellar populations significantly depends on the flux den- 
sities reached by the radio survey. In the lower panels of fig- 
ure |9] the sample plotted in the upper panel has been split 
in order to draw more meaningful conclusions on the stel- 
lar populations of flux-limited radio-selected sub-samples. 
As it is evident from figure |9l and as it is expected (see sec- 
tion [T]), the stellar population properties of the host galaxies, 
and thus likely also the process responsible of the radio emis- 
sion, are different in different radio flux density ranges. In the 
shallowest sample considered {Siachz > 84/iJy), the fractions 
of blue, green, and red classified objects are roughly similar 
(33%, 36%, 30%, respectively); in the Siaghz > 40/iJy and 
SiAGHz > 24/iJy samples they are about 40%, 40%, 20% , and 
eventually they become 45%, 33%, 22% in the deepest ra- 
dio sub-sample (Siaghz > 16/iJy). This clearly suggests how 
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Fig. 10. — The age/r distribution for the sub-samples of quiescent, inter- 
mediate, and star-forming objects. The solid symbols and errorbars on top 
of the histograms show the median and 16'''-84''' percentiles of each distribu- 
tion. Age/r values larger than 9.5 were set to 9.5 for the purpose of this plot. 
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Fig. 1 1 . — The relative contributions of different SED-clas.sified sample.s 
as a function of radio flux density. The points show the median 20cm flux 
density in the flux bin (indicated by the horizontal error bar) against the frac- 
tion over the whole population in that flux bin. Color coding is the same as in 
figure [9] The flux density bins for the different SED-classified samples are 
the same, but are shown with a small offset in the figure for clarity. 

the contributions of the actively star-forming and red galaxies 
change with the limiting flux density of the sample, with red 
galaxies increasing their relevance in brighter samples and, 
viceversa, star-forming galaxies becoming more important in 
fainter samples. For these flux-limited samples, according to 
a Kolmogorov-Smirnov test, the U-B color distributions of 
the radio-faintest and brightest sub-samples are different at 
a 2± 98% significance level. However, the significance of the 
change in the colors of the host galaxies at different radio flux 
densities is obviously more evident when considering sources 
in ranges of radio flux density instead of flux-limited sam- 
ples, as we will do in the following. 

We should notice that, since the lower panels of figure 
|9] refer to a flux-limited sample, there are different effects 
which come into play in realizing the different distributions 
of the different-depth sub-samples, as for instance the higher 
fraction of z ~ 1 as compared to z < 0.5 objects probed 
at fainter flux densities (see figure |6]l, and likely an evolu- 
tion with redshift of the 20cm luminosity threshold between 
star formation-dominated and AGN-dominated radio sam- 
ples (see figure|8]l. This is further discussed below. 

The different (sky-coverage corrected) contributions of the 
red, green, and blue classified sources are shown in figure [TTI 
This figure shows how the S\,4ghz ^ 100/iJy sources, which 
have a relatively flat redshift distribution in the redshift range 
we are probing (fig. |6]l, are roughly equally split between the 
three sub-classes of red, green and blue galaxies. When go- 
ing to fainter flux densities, where we saw that the redshift 
distribution becomes more skewed toward higher redshifts, 
figure [TT] shows that the sample is depleted of red galax- 
ies while the fraction of blue actively star-forming systems 
increases (green-classified objects make up ss30-40% in all 
sub-samples). In the faintest subsample, with Siachz in the 
range 16-25/iJy, less than 25% of the 0.3 < z < 1.3 sample is 
at 0.3 < z < 0.7, compared to more than 40% in the shallow 
SiACHz ^ 100/iJy subsample; the faintest sample only hosts 
^10% red galaxies, while the fraction of blue galaxies has 
increased to more than 50%. According to a Kolmogorov- 
Smirnov test for these radio-faintest and brightest samples, 
the U-B color distributions are different at a >99.7% signifi- 
cance level. 



We also note that, by splitting host galaxies just based on 
their U-B restframe color (e.g., red galaxies with U-B>1, 
blue galaxies with U-B<1) we find that , in agreement with 
previous work (e.g. iMainieri et al.l l2008h . the radio popula- 
tion is dominated by red galaxies above flux densities of 
100/iJy, while below 80/iJy blue galaxies begin to domi- 
nate (« 65%, considering a flux-limited sample with 40 < 
S\ 4CH- < 80/iJy). Splitting this sample in three redshift bins 
with 6.3 < z < 0.6, 0.6 < z < 0.9 and 0.9 < z < 1.3, we find 
that this « 65% fraction stays basically constant as a func- 
tion of redshift. However, we should think in terms of lu- 
minosities instead of flux densities: while the flux density 
range 40 < Siachz < 80/iJy at z ~ 0.3 corresponds to a lu- 
minosity range of approximately 1-3 x lO^^W/Hz, at z 1.3 
it corresponds to luminosities of order 3-6 x lO^^W/Hz. These 
^lO^^W/Hz luminosities in the lowest redshift bin would cor- 
respond to flux densities well above lOO/xJy, which as we said 
are usually dominated by red galaxies. Therefore, the appar- 
ently mild evolution between redshift z ^ 0.5 and z ^ 1.1 of 
the fraction of blue galaxies in Siachz < 80/iJy samples ac- 
tually happens against a change in radio luminosities of an 
order of magnitude. This is due to evolution in these faint 
populations as we will further discuss below ( section |7!2] |. 

In the following we compare our SED-selected subsamples 
in terms of different classification criteria from previous stud- 
ies. 

7.1.1. NUV/optical vs radio properties 

For the sake of comparison, we show in figure [12] how our 
classification based on the observed broad-band SED com- 
pares to two different classifications of radio samples based 
on specific rest-frame properties. In the top panel of fig- 
ure [12] we_com£^ewith_flie AGN/star-forming classification 
used in lSmolcid et"an d2008h . Based on a SD SS/NVSS/IRAS 
sample of local radio sources, ISmolcic et all (2006 2008) de- 
vised a classification method to separate populations of ra- 
dio sources whose 1 .4GHz emission is dominated by AGN, 
by star formation, or is likely to be contributed by both pro- 
cesses. This is based on the principal component restframe 
colors PI and P2, which are linear combinations of restframe 
colors in the modified Stromgren syste m in the w avelength 
ran ge 3500 - 58 00 (Odell et al. 2002, Smolcic et al. 2006). 
In Sm olcic et al.l (^2008), a color cut at Pl=0.15 was adopted 
to separate their radio sample in the COSMOS field into two 
populations of AGN and star-formation dominated systems. 
As the figure shows, indeed also our objects with a PI > 0. 15 
are mainly classified as quiescent (thus likely with an AGN 
produced radio emission), or at most as intermediate color 
sources (thus with a possibly relevant contribution by AGN). 
On the other hand, our class of intermediate objects extends 
down to slightly below P1=0, where increasingly star forming 
populations become dominant extending down to Pl« —0.7. 
Thus, to first approximation, most of our green objects would 
be classified as star- forming in this scheme^^. Nonetheless, 
from e.g. Smolcic et al. (2006), we note that a significant frac- 
tion of the sources with PI around zero (roughly —0.1 <P1< 

We note thatlSmolcic et al.' C008f) applied a coirection to the synthetic 
PI color used for the COSMOS sample, through comparison of the PI colors 
estimated by the SED-fit model and by the spectrum for a sample of SDSS 
sources. No such connection was applied here. If applying to our PI colors the 
same corrections used in Smolcic et al. (2008), the PI color of green objects 
would be on average lower by about 0.06 mag, thus moving further into the 
range of sources classified as star-forming. 



12 



Strazzullo et al. 





11 12 13 14 15 

Log(L,.,cH. / M.) 

Fig. 12. — Top panel: The principal restframe color-color diagram (P1,P2) 
for objects with redshift 0.3 < z < 1.3. Symbols and color coding are based 
on the best-fit SED type as in figure|2] The solid black line shows the color- 
cut Pl=0.15 b etween AGN and star-forming dominated populations used in 
ISmolcic et all (2008). Bottom panel: The A. (4000) index plotted against 
the radio luminosity per unit stellar mass, for objects with redshift 0.3 < z < 
1.3. Symbols and color coding as in figure [9] The dashed line s hows the 
separation between AGN and stai-forming galaxies as defined in Best et 
j2005b). For the purpose of this plot stellar masses were scaled to a Krou 
BOOK) IMF. 



0.1) populates the region in the BPT ([Baldwin et al.ll 198 1*) di- 
agram where the AGN/star-forming classification is consid- 
ered uncertain. 

In the bottom panel of figure [12] we plot instea d our sam- 
ple in t he LogiLi^cHz 

/M*) vs A, (4000) plane. iBesF et alj 
( I2005bl) divided a sample of SDSS/NVSS matched sources 
at redshift z < 0.3 into two broad classes of AGN and star- 
bursting galaxies. The selection was in fact based on the 4000 
break index D„(4000) and radio luminosity per unit stellar 
mass LiACHz/M^,. Since we don't have spectra for all sources 
in our sample, we use synthetic D„(4000) indices derived 
from best-fit stellar population synthesis models, which may 
only be considered as a very approximate estimate of the real 
D„(4000). We use the stellar masses described above to cal- 
culate Li,4GHz/M* for each galaxy. If the estimates we plot are 
representative enough of the true D„(4000) and Lj 
we should conclude that indeed basi cally all of our blue ob- 
jects are below the Best et alJ (12005 b) division line, and the 
vast majority (>90%) of the red galaxies are above the line^**. 

As far as these minor differences are concerned, we note that the exact 



The green sources fill the gap between the two, just about the 
division line, falling in both the "AGN" and "starburst" re- 
gions. 

Finally, we show in figure [13] how our classification of this 
sample compares with other galaxy populations with regard 
to radio/UV flux densities. Figure [13] shows for our sam- 
ple the star formation rate (SFR), as determined from the 
1 .4GHz luminosity, against the dust attenuation estimated as 
Ajgoo =2.5Log(SFR,.fl,/,o/SFRf/y). We note that the meaning 
of these quantities for the whole sample is ill defined. In fact, 
only for galaxies whose radio emission is due to star forma- 
tion the quantity plotted as "SFR" actually represents the SFR, 
and A|g()Q indeed is the dust attenuation. Instead, galaxies 
hosting a radio emitting AGN have their 1.4GHz luminos- 
ity at least partially contributed by the AGN, and thus both 
"SFR" and Aj^qq lose their meaning for these objects. How- 
ever, keeping this in mind, this plot allows us to compare in 
a simple way our sample with other populations of star form- 
ing galaxies at different redshifts and selected with different 
criteria^^. 

As the figure shows, the range of attenuations derived for 
our blue objects is in very good agreement with other stud- 
ies of different kinds of star forming galaxies at different red- 
shifts iC alzetti et al. 2000, Calzetti 2001, Hopkins et al. 200 ll 
lAfonso et al.. ,2003. .Choiet al-: .2006., Pannella et al.. .2009I) . 
While part of the green objects would also overlap with these 
samples, it is clear how the Ajgoo for the green population 
generally lies above the expectations, and definitely the red 
population has too high values of A*goQ. This might support 
the idea that the red population is for the great majority made 
of AGN hosts, and also suggest that at least part of the green 
population has a contribution from AGN to the radio luminos- 
ity, even though part of these galaxies can still be very dusty 
systems. 

The significant occurrence of AGN hosts among the in- 
termediate population between red and blue g alaxies has 
been noted in several previous studies ( e.g., IChoi et aT] 
l2009l iMartinltallllOOTl [Nandra et al.ll200l and references 
therein), as well as their actual nature of composite (meaning 
SF-hAGN) system s (jSchawinski et al. 2007, Wild et al. 200"^ 
iSalim etaW flMf). We note that the classification of such 
composite sources in the literature is quite variable, and for 
instance galaxies with more than 10% of their radio luminos- 
ity contributed by s tar form ation have been classified, in some 
cases, as starbursts dTasse et al. 2008). This may indeed be the 
case, and certainly also in our green sample different amounts 
of star formation are presents, however a starforming-or- 
AGN classification may be not appropriate for these sources, 
especially depending on the kind of study they are used for 

7.1.2. Radio— IR properties 

As discussed above, in spite of the insight that we can cer- 
tainly gain by combining radio fluxes and optical/NlR pho- 

shape of the division line was determined, also based on emission-line diag- 
nostics, for a local (z < 0.1) galaxy sample, while the galaxies in our sample 
were observed 2 to 7 billion years earlier, as it is clear from the D„ (4000) 
range. 

We note that the plotted P^2iOQ ^'^'^ values from other studies have 
been derived from the original published quantities with some assumptions, 
namely: 1) the corrected SFRs derived from 1.4GHz, Hq, UV or IR lumi- 
nosities agree with each other; 2) the dust obscuration estimated from the 
comparison of radio or IR vs UV fluxes, UV slope, or Balmer decrements, 
also agree with each other once the appropriate translations are made; 3) the 
Calzetti et al. (2000) law; and 4) the color excess for the stellar continuum is 
a factor 0.44 of the color excess for the nebular gas emission lines. 
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Fig. 13. — The "dust attenuation" A*gQ(,=2.5Log(SFR„rf,o/SFR28oo) 
against "star formation rate" as determined from tlie 1.4GHz luminosity 
(see text for tlie actual meaning of the two quantities plotted). Small sym- 
bols color coded as in figure |9] show the 0.3 < z < 13 radio sample (21 
objects (out of which >95% are red or green) with A|gQQ > 7.5 and/or 
SFRrudio > 2000 are not shown). Large empty circles with errorbars show the 
IPan nella et al. (2009) results (median and 16-84''' percentile range) for star 
forming BzK galaxies atz~2. Empty squares and stars show local galaxies 
from Ca lzetti et aU (2000). Calzetti (2001) (star symbols are for galaxies clas- 
sified as dust-^rich in Calzetti (2001); note NCG 1614 with SFR=55M0/yr 
and very high dust extinction, which by the way hosts an obscured AGN (e.g., 
[Guainazzi & Bianchi (2007))). The solid, dotted and dashed lines show the 
dust attenuation as a function of SFR as derived from the local optical/UV- 
selected sample of Hopkins et al. ( 2001), the intermediate redshift NIR/MIR- 
selected sample of Ch oi et al. ( 2006) and the low-redshift radio-selected 
sample of lAfonso et alj i2003ll . respectively (see text for details). 

tometry, the results presented so far may not provide conclu- 
sive evidence about the actual nature of these sources. There- 
fore, we need to introduce further information which may help 
us identify which process powers the radio emission. Obvi- 
ous promising data already available on this field are X-ray 
(Chandra) and infrared (Spitzer and Herschel) observations. 
While we postpone a full analysis of these data to a future 
work, we use here just the Spitzer/MlPS 24/im data to es- 
timate the total infrared (IR) luminosity of our sources and 
thus examine the behavior of our SED-selecte d sub-samples 
with respect to the FlR-radio correlation (e.g.. lCondonlll992l 
lYun et alj|200lb . The DSP was observed at 24^m with MIPS 
onboard Spitzer as part of the GO-3 program #30391 (PI: 
F. Owen), for a total of 60.6 hours over an area of about 
half square degree, and a median integration time per pixel 
of about 2500s. The 5a flux density limit is estimated to 
be about 40/iJy. The data reduction and catalog production 
are described in full detail in a forthcoming companion paper 
(Owen et al., in preparation). More than 80% of the faint radio 
sources from the 90% complete sample in the redshift range 
0.3 < z < 1 .3 are detected at 24/im. Matched sources with flux 
possibly contaminated by neighbors within the 24/im PRE are 
estimated to be about 10% of all matched detections, based 
on the cross-correlation with the IRAC 3.6/im catalog. While 
upper limits make up for < 20% of the whole 0.3 < z < 1.3 
sample, they are more relevant for red-classified objects (al- 
most 40% of upper limits) than for green and blue sources. 

In the redshift range 0.3 < z < 1.3, the observed 24/tm light 
probes the restframe ^ 10 — 20/im. We use the templates of 
ICharv & ElbazI ( 1200 Ih to estimate from the 24/im flux density 
the total (8-1000/im) IR rest-frame luminosity (or an upper 
limit for 24//m-undetected sources). This is done taking the 



template whose predicted luminosity at the observed 24fim is 
closest to the actual observed 24/im luminosity. Taking the 
mean or median luminosity across the whole template library 
typically changes this number within about 0. 15dex, but since 
this is obviously library-dependent we don't see a real point 
of adopting the mean or median instead of the formal best- 
fitting template. Furthermore, even just based on the tem- 
plates available in this library, such estimate of the total IR 
luminosity based on the single observed 24/im point may be 
affected by systematics of up to a factor 2. 

We use this total IR luminosity together with the 1.4GHz 
luminosity to estimate the logarithmic ratio of bolometric IR 
and monochromatic radio luminosities ^/R=Log(L/fi/(3.75 • 
10'^)) - LogiLiACHz) (lHelouetal.lll985h . This is plotted as 
a function of redshift for the different sub-samples of red, 
green, and blue sources, in the bottm-left panel of figure 
[l4l In this figure, filled circles show sources unambigu- 
ously matched with a 24/im source, and filled triangles show 
sources matched within 2" with a 24/im source but whose 
24/tm flux might be contaminated neighbors. Finally, upper 
limits are shown by down-pointing arrows. 
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Fig. 14. — Bottom panels: in the left panel the flux ratio qiR is plotted 
against redshift for the different SED-selected sul3-samples (color coding as 
in figure|9). Filled circles show unambiguously matched 24/i sources, filled 
triangles show matched 24/i sources whose flux might be contaminated by 
neighbors, and arrows show upper limits (see text for details). Large orange 
symbols show m atched X-ray s ources from the Trouille et al. (2008) sample 
(according to the lTrouille et alj {2008) spectral classification, circles, down- 
ward triangles, triangles and squares show absorbers, star formers, high- 
excitation sources and broad-line AGNs, respectively). The gray-shaded 
area shows the la range about the (redshift independent) tjiR = 2.4 reported 
in Ivison et al. (2009). The right-hand panel shows the sky-coverage cor- 
rected distributions of y/fi for the sub-samples plotted in the left panel, in- 
cluding upper limits (see text for details). The histograms for the red, green, 
and blue samples are evaluated in the same qiR bins, but are shown slightly 
offset for clarity. Top panel: same as bottom left panel, but for the oh- 
served (non k-corrected) flux density ratio ^24- T he gray sha ded area shows 
the envelope of Chary & Elbaz ( 200 1 ) and Rieke e t alj 120091) templates with 
L//i=lO"L0 and W^-Lq (plus IO'^Lq for z < 0.5). 

The bottom-right panel of this figure shows the distri- 
bution (corrected for the sky-coverage) for the three different 
red, green, and blue sub-samples. It is important to note that 
the plotted histograms (as well as the related statistics given 
below) include both 24/tm detections and upper limits. This 
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does not affect our results (or actually affects our conclusions 
in a conservative way), since we are interested in the differ- 
ence between the qiR distribution of the three sub-samples, 
with qiR of intermediate and quiescent sources expected to be 
lower than of star forming sources, if the radio power of 
objects classified as intermediate and quiescent is at least par- 
tially provided by an AGN. This means that, including upper 
limits (and possibly contaminated 24/im detections), we are 
- if anything - reducing the actual difference between the qiR 
distributions of the three sub-samples. 

In both panels, the gray-shaded area shows the Icr 
range about the (re dshift independent) <7/« = 2.4 reported in 
llvison et all ( l2009l) . As fi gure [14] shows, blue star-forming 
sources mostly he on the expected FlR-radio correlation. 
Also up to 50% of the green sources lie within Icr of the 
expected FlR-radio correlation, which would point toward 
these being powered by star formation as well. However, 
even though with a broader distribution, red sources also lie 
close to the FlR-radio correlation, with up to 40% of this 
sub-sample lying within Ict of the correlation. Sources pop- 
ulating the gray-shaded area in figure [14] ( 1 cr about the ex- 
pected correlation) are for more than 50% classified as star- 
forming, but show a sizable fraction of ^ 30% and ^ 10% 
of green and red sources, respectively. The distribution 
of red sources appears to be very different from those of blue 
and green sources. Even though a Kolmogorov-Smirnov test 
cannot be meaningfully applied to the sky-coverage corrected 
distributions, applying it to flux-limited sub-samples selected 
in portions of the radio image with depth uniformly better 
than a given threshold, suggests than the qiR distribution of 
red sources is different at a significance of more than 99.9%. 
A test on the binned distributions (including errors) plot- 
ted in figure[T4](as well as on similar distributions binned with 
half bin size), also suggests that the distributions of qiR of red 
vs. blue or green sources are different at a > 99.5% level. On 
the other hand, the qiR distributions of blue and green sources 
look much more similar The test on the binned distribu- 
tions suggest that they are different at a > 98% level, and the 
Kolmogorov-Smirnov test on a limited part of the sample (as 
above) gives a Pks ^ 0.2, thus based on the present data these 
distributions are, at most, marginally different. 

Based on these results one might thus conclude that a large 
fraction of our sources with intermediate colors, as well as 
a sizable fraction of the red ones, are actually reddened star- 
bursts. On the other hand, we also note that among the sources 
which are X- r ay det ected (and mostly AGN classified in the 
iTrouille et al.l (l2008h sample), many lie on the FlR-radio cor- 
relation, including some classified as absorbers (based on Oil 
or Hq,h-NII equivalent widths). Indeed, other studies have 
found that sources classified as low-radio power AGN may 
often lie on or close to the same FlR-radio correlation ex- 
pected f or star formin g sources (e.g. see discussions and refer- 
ences in lSmolcic et al . 2008, Sargent et al. 2010). However, it 
should also be noted that an AGN selected based on its (spec- 
tral or photometric) optical/NlR properties (or X-ray), hosted 
in a "composite" system with ongoing star formation, may be 
in a stage where radio emission is negligible and thus does 
not significantly affect the IR/radio properties which remain 
determined by the star formation process. In this case, while 
information at other wavelengths suggest the presence of an 
AGN, it might not (significantly) contribute to the radio power 
of the source. Furthermore, obscured star formation confined 
in limited areas of a galaxy might also determine its IR/radio 
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Fig. 15. — The 1.4GHz luminosity function of different galaxy classes. 
Large empty symbols refer to the whole populations, while filled symbols 
and their coding refer to the same classes as in figure |9] Number densities 
are shown in three redshift ranges, as indicated. The luminosity functions for 
the different samples are evaluated in the same L\ ^Q^j- bins, but are shown 
slightly offset for clarity. Errorbars show Poissonian errors. 

properties while possibly going undetected in broad-band 
photometry at optical wavelengths. On the other hand, we re- 
mind again the reader that the qiR in figure [T4]relies uniquely 
on the observed 24pim flux density, and thus might be biased 
producing a spurious result. For comparison, we also show 
in the top panel of figure [14] the observed (non k-corrected) 
flux density ratio ^24 — Log{S24fim/S\.4GHz)- As a reference 
the gr ay s haded area shows the envelope of IChary & ElbazI 
(120011) and lRieke et"aD(l2009h templates withL/R=10"Lo and 
IQ^^Lq (plus also IQ^^Lq for z < 0.5). This ^24-based figure 
essentially reflects the results of the bottom panel. We will 
continue the investigation of IR/radio properties of this sam- 
ple with a proper SED analysis of the full data in a future 
work. 

7.2. Luminosity functions and evolution of faint radio 
populations 

Finally, we investigate the contributions of the different 
galaxy classes to the sub-mJy population by plotting in fig- 
ure \T5\ the 1 .4GHz luminosity functions (total and split by 
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Fig. 16. — The redshift evolution of the 1.4GHz luminosity function of 
galaxies in different SED classes, as indicated in each panel. Luminosity 
functions are the same as those plotted in figure [Ts] Different colors refer to 
different redshift ranges, as also indicated. Errorbars show Poissonian eiTors. 
Small squares (with same color coding as so lid lines) show, as a reference, 
the LF of radio-selected AGN as derived bv lSmolcic et alj 42009bl) in similar 
redshift bins (see text for details). 

SED class) in three redshift ranges. Figure [15] shows num- 
ber densities based on the 0.3 < z < 1.3, 90% complete 
sample, estimated wi th the l/Vma.x method ( ISchmidtlll968l 
lAvni& BahcaIilll980h . Within the redshift range 0.3 < z < 
1.3, the maximum volume over which an object can enter 
our sample essentially depends on its radio luminosity and 
on the (non uniform) depth of the 1 .4GHz image. Therefore, 
for each source of a given radio luminosity Li,4g//;, the max- 
imum volume accessible to the source was calculated, based 
on the maximum redshift out to which the source would have 
been detected as a function of the varying image depth of the 
1.4GHz image. We remind the reader that, similarly to what 
discussed above concerning the sky coverage correction, the 
l/Wmnx correction is not negligible for low-luminosity sources 
due to the limited survey area probing the faintest fluxes. On 
the other hand, we also note here that luminosity functions 
obtained, without the use of the IN max correction, by using 
smaller volume limited samples defined in portions of the ra- 
dio image and luminosity ranges, are generally in very good 
agreement with those determined with the l/Vmnx method, and 
lead to the same conclusions. 

Only sources in the > 90% complete sample were used, 
and no further correction was adopted for the small residual 
incompleteness of the sample. However, we remind the reader 
again that the > 90% completeness level of this sample was 
estimated based on some assumptions (as discussed in section 
5.3). 

Figure [15] shows the relative contributions of the differ- 
ent SED-classified sub-samples as a function of redshift and 
luminosity. It also suggests an evolution of the luminosity 
functions of all classes over the probed redshift range, which 
is shown more clearly in figure [16] This is confirmed by 
a fit to the binned data in the thr e e red shift bins with 
a parametric LF of the ISaunders et al.l (11990 ) form $ (L) = 
$*(L/L*)i-"exp(-0.5(T-2[log(l +L/U)]^). We note that, 
in order to keep consistency within our SED classification, the 
fit was performed only in the three redshift bins shown in fig- 
ures[T5]and[T6] and over the luminosity range properly probed 
by our data (as a reference, Log(Li 4chz)< 24.1,24.4,24.6 in 
the z ~ 0.5,0.8,1.1 redshift bins, respectively), without in- 



cluding measurements from other surveys sampling brighter 
luminosities, nor a z=0 reference LF. Nonetheless, we show 
in figure [16] the LF of radio selected AGNs as derived by 
[Smolcic et al. (2009b) in the redshift bins 0. 35 <z< 0.6 
0.6 < z < 0.9, 0.9 < z < 1.3. We note that the 'Smol5ic et alj 
( l200 9b) LFs are shown just as a reference, as the AGN sam- 
ple of Smolcic et al. does not perfectly match our red sam- 
ple, as shown in figure [12] top paneF . Given this differ- 
ence in the sample selection (and a small difference in the 
first redshift bin), the small area of our field which does not 
properly probe the brighter luminosities better sampled by the 
large COSMOS survey, and in general the uncertainties in- 
volved in the LF determination, our red-sample LFs can be 
co nsidered in rea sonably good agreement with the AGN LFs 
of lSmolcic etai] (l2009b). 

For each SED class (red, green and blue), as well as for 
the total population, we performed a simultaneous fit to the 
binned LFs in all three redshift bins, allowing for redshift evo- 
lution in the form L oc (Ih-z)"'^. This assumption of pure lu- 
minosity evolution (PLE) is only used as a working tool in 
order to quantify the observed evolution, and for comparison 
with other studies. In fact, we are well aware of the fact that, 
beside the well known degeneracy between luminosity and 
density evolution, there is indeed little reason to believe that 
either pure luminosity or pure density evolution (PDE) may 
be an adequate descri ption over a range of different luminosi- 
ties and redshifts (e.g jDunlop & Peacoc5ll990irMllott et alj 
HOOI, Ueda et al. 20031) • The very fact that radio populations 
are made of different types of objects, and that these different 
sub-populations might evolve in different ways, imply that 
simple PLE or PDE cannot be in general considered as an 
adequate description. On the other hand, due to the limited 
luminosity range probed, and to the small number statistics, 
our data alone cannot be sufficient to effectively constrain the 
general 20cm luminosity function and its redshift evolution, 
so we will only try to quantify the amount of evolution ob- 
served, at the luminosities and redshifts properly probed by 
our data, by assuming the simple PLE model of redshift evo- 
lution. 

The best-fit ai obtained for the whole population is 3.5 ± 
0.2, while for the red, green, and blue sub-samples we ob- 
tain a/,=2.7±0.3, 3.7+„4, and 3.2+„2' respectively. A non- 
parametric evaluation of a^, obtained by directly comparing 
the LFs in the three redshift bins (and assuming L ex. (1-Hz)"'- 
as above), instead of fitting the parametric Sau nders et alj 
(1990) form, yields results perfectly consistent with those 
listed above for the parametric fitting (3.1 ±0.2, 2.5 ±0.3, 
3.6 ±0.2, 2.9 ±0.3, for the total, red, green, and blue sam- 
ples, respectively). 

The formal best-fit ai is thus close to w 3 for the whole 
sample as well as for all SED sub-samples. This suggests 
that our observations are consistent with luminosities decreas- 
ing by a factor « 10 from redshift just above 1 to the local 
Universe. The fact that the evolution factors for the three 
SED classes are very close to each other might suggest also 
in this work a link between the evolution of star formation 
and AGN activity (e.g., Silverman et al. 2008, Heckman 20()^ 
and references therein), provided that our SED-selected sub- 
samples actually probe such different populations. 

A PLE rate az. « 3 is similar to PLE rates esti- 

Because of the even greater difference between our blue and green sam- 
ples and the Smolcic et al. star-forming sample we do not attempt a compar- 
ison of the LFs for star-forming galaxies. 
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mated in other studi es for star-forming galaxies (e.g., 2.5 
dSevmour et al.ll2004). 3 j Cowie et al. 2004), 2.7 with a negli - 
gible an=0.15± 0.6 (Hopkins 2004), 2.7 (Huvnh et al.!l2005|), 
2.3 (iMoss et al. l 2007), 2.1-2.5 (Smolcic et al. 2009a)) or X- 
ray AGNs (e.g., 3.2 (Barger et al. 2005), 2.7 (Delia Ce ca et alJ 
HOOS)). As far as AGNs are concerned, we should note though 
that low-luminosity radio AGNs have been found to show 
slower evolu ti on com pared to higher radio power AGNs (e.g., 
IWillott et alj (120011) ) and, at luminosities similar to those 
probed here, somewhat lower PLE rates as compared to our 
red sample have been measured in some previous work (e.g. 
aL=2 ( Sadler et al. 2007) or a£=0.8 (ISmolcic et all2009bi) us- 
ing the .Sadler et al.. (2002) AGN LF). 

8. SUMMARY 

We have carried out a multi-wavelength analysis of faint 
radio sources in the Deep Swire VLA Field. The depth of 
this survey allows us to probe populations of radio sources 
uncommon or absent in other deep radio surveys, with almost 
a thousand sources fainter than 50/iJy. 

Based on optic al/NlR/lR AC photometry, we built the SEDs 
of the identified counterparts and compared them both with a 
galaxy SED template library and with stellar population syn- 
thesis models, determining their photometric redshifts, stellar 
masses, and broad stellar population properties. The derived 
redshift distribution of radio sources appears to be different 
at different flux density levels, with the distribution of the 
faintest sources showing a more pronounced peak, at about 
redshift one. 

We have focused on a 90% complete sample of counterparts 
of sub-mJy radio sources with redshift 0.3 < z < 1.3, divid- 
ing the sample in broad classes of quiescent, intermediate and 
star-forming systems based on their optical/NlR colors. The 
population mix as described by these sub-samples shows a 
clear dependence on radio flux density, with an increase of 
star-forming populations at lower fluxes, in agreement with 
previous studies. At all redshifts up to z 1.5, the contribu- 
tion of star-forming galaxies becomes increasingly important 
at lower radio luminosities. 

The rest-frame U-B vs B color-magnitude diagram of this 
radio-selected sample shows the presence of a significant 
green-valley population, beside two populations of red and 
blue galaxies. One might assume that the radio emission from 
red galaxies is in most part due to an AGN, because of their 
apparently very low star formation (beside a possible contri- 
bution of extremely dusty galaxies), and on the other hand 
that there is a predominant contribution to the radio emission 
from star formation in blue star forming galaxies. 

The properties of the intermediate, green valley sample 
with respect to the "AGN" vs "star-forming" classification 
are less clear. Their stellar population properties, and the 
comparison of radio and UV luminosities, suggest that this 
class of objects might be a mixed population of AGN and 
star-forming galaxies, possibly including composite systems 
where both nuclear activity and star formation are present. 
On the other hand, the comparison of 24/im-based IR and 
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1.4GHz radio luminosities would suggest that many of our 
faint radio sources, and not limited to the blue, obviously 
star-forming galaxies, may lie close or on top of the radio- 
IR correlation expected for star forming objects. In particular, 
not only the radio-to-lR flux ratios of sources classified as 
intermediate are distributed very similarly to those of star- 
forming galaxies, but also up to 40% of the "quiescent" clas- 
sified galaxies have a measured in the range 2.2-2.6. This 
result may be affected by the fact that we are sampling the IR 
SED only with one (24/im) photometric point. More stringent 
conclusions will be possible thanks to the analysis of the full 
IR SEDs, which will be the subject of a forthcoming work. 

In agreement with previous studies, all populations studied 
in this work show evolution with redshift. In a simple PLE 
scenario, 1.4GHz luminosities decrease by about a factor 10 
from redshift just beyond one to the local Universe. 

A comparison with the whole sample of radio-undetected 
objects in the field, as well as a stacking analysis to study av- 
erage radio properties of complete mass/SED-selected galaxy 
samples, are postponed to a future work, and will likely lead 
to a better understanding of the nature of the populations stud- 
ied here, as well as more in general of the interplay of AGN 
and star formation activity in the evolution of galaxies. 
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TABLE 1 

Values of the magnitudes maGc, and magiq derived through simulations of photometric measurements in apertures of diameter 
1 .5" in dual-image mode of artificial point-like sources added to the images. the images marked with a * were smoothed to a 

FWHM=1.3". See text for details. 



passband 


magn,, 


magio 




(AB mag) 


(AB mag) 


NUV 


25.5 


24.0 


U 


28.0 


25.5 


g* 


27.8 


25.3 


r* 


27.3 


24.9 


i* 


26.5 


24.3 


z* 


26.2 


23.6 


J* 


25.2 


23.0 


H* 


25.0 


22.7 


K* 


24.2 


22.4 


3.6/im 


24.4 


22.4 


4.5fim 


24.2 


22.2 



TABLE 2 

Some of the main derived properties used in this work. This table is available in its entirety on the electronic edition, a portion 

IS SHOWN here for GUIDANCE. ONLY OBJECTS WITH AVAILABLE SPECTROSCOPIC REDSHIFT OR RELIABLE PHOTO-Z ARE LISTED. SEE TEXT FOR 
DETAILS ABOUT HOW THESE QUANTITIES WERE DERIVED. NOTES: (") ID, RA AND DEC AS IN PAPER L (*) SPECTROSCOPIC REDSHIFT WHERE 
AVAILABLE, OTHERWISE PHOTOMETRIC REDSHIFT. REDSHIFTS WITHOUT AN ERROR LISTED ARE SPECTROSCOPIC. C) QUALITY FLAG FOR 
PHOTOMETRIC REDSHIFTS, AS DEFINED IN SECTIOn ISTTI ('') SYNTHETIC ESTIMATES OF ABSOLUTE MAGNITUDE IN THE R (JOHNSON) BAND, 
RESTFRAME U-B COLOR, AND STELLAR MASS (SALPETER IMF), ASSUMING THE REDSHIFT LISTED IN COLUMN 4. BOTH R BAND MAGNITUDE AND 
STELLAR MASS ESTIMATES ARE BASED ON SED FITTING IN A I .5"APERTURE AND CORRECTED FOR EACH OBJECT BASED ON THE DIFFERENCE 
BETWEEN TOTAL AND APERTURE MAGNITUDES AS MEASURED IN THE DETECTION IMAGE, AS EXPLAINED IN SECTI0n[6] 



ID" 


RA" 


Dec" 


redshift* 


^^photo — z 


K 


U-B'' 


Log(M*)'' 




(J2000) 


(J2000) 




(mag) 


(mag) 


(Mo) 


20032 


10:43:41.58 


59:09:20.4 


1.66 ±0.50 


AA 


-20.9 


0.55 


10.35 ± 0.25 


31 


10:43:43.59 


59:00:10.2 


0.92 ± 0.05 


AA 


-21.9 


0.84 


10.26 ± 0.20 


35 


10:43:43.93 


58:56:51.6 


0.64 ± 0.07 


AA 


-21.3 


0.79 


10.09 ± 0.26 


20038 


10:43:44.33 


59:05:51.1 


0.96 ±0.10 


AA 


-21.1 


1.01 


9.96 ± 0.32 


36 


10:43:44.76 


59:15:02.9 


1.62 ± 0.09 


AA 


-24.2 


1.09 


11.00 ±0.22 


37 


10:43:44.84 


58:58:18.7 


0.52 ±0.11 


AA 


-22.8 


1.17 


11.24 ±0.13 


10056 


10:43:45.22 


58:57:40.2 


2.08 ± 0.24 


AA 


-22.6 


0.39 


10.26 ±0.21 


20045 


10:43:46.13 


59:01:37.9 


1.04 ±0.07 


AA 


-21.6 


0.93 


10.24 ± 0.39 


44 


10:43:46.21 


59:01:19.1 


1.34 ±0.08 


AA 


-23.9 


1.01 


11.02 ±0.14 


46 


10:43:46.91 


59:00:21.3 


0.98 ± 0.05 


AA 


-21.9 


0.98 


10.35 ± 0.36 


10061 


10:43:47.91 


59:06:21.2 


0.14 




-19.3 


0.57 


9.60 ± 0.26 


47 


10:43:48.05 


59:02:23.2 


2.26 ± 0.45 


B 


-23.8 


0.75 


10.90 ± 0.22 


49 


10:43:49.22 


58:55:38.7 


1.58 ±0.33 


AA 


-21.0 


0.44 


9.59 ± 0.28 


20053 


10:43:49.26 


58:51:34.7 


0.74 ±0.10 


AA 


-22.0 


1.10 


11.03 ±0.18 


50 


10:43:49.46 


58:56:31.5 


1.06 ±0.10 


AA 


-22.5 


0.97 


10.69 ±0.18 



